Gravitational waves from inspiraling compact binaries: 
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We present further evidence that the second post-Newtonian (pN) approximation to the gravita- 
tional waves emitted by inspiraling compact binaries is sufficient for the detection of these systems. 
This is established by comparing the 2-pN wave forms to signals calculated from black hole pertur- 
bation theory. Results are presented for different detector noise curves. We also discuss the validity 
of this type of analysis. 
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I. INTRODUCTION 

There is currently a large effort underway to build 
an international network of interferometric gravitational 
wave detectors. These include the two American LIGO 
detectors, the French-Italian VIRGO detector, the some- 
what smaller British-German GEO600 and the Japanese 
TAMA detector. 

One of the most promising sources for detection are 
compact inspiraling binary systems |Q]. These systems, 
consisting of a pair of black holes or neutron stars, lose 
energy through the emission of gravitational radiation 
and therefore spiral inward. The emission is strongest 
during the last few minutes of the inspiral; just before 
the orbit becomes unstable and the two objects plunge 
towards each other. For stars with masses a few times 
that of the sun this last burst of waves falls within the 
frequency band of the previously mentioned detectors. 

The method of choice to search for such signals in a 
noisy data stream is matched filtering g] , where the de- 
tector output is compared to a bank of templates. A 
good overlap or match then indicates the presence of a 
signal. It can be shown Q that matched filtering is an 
optimal method. However, to use this technique an accu- 
rate knowledge of the expected signals is necessary to set 
up the template bank. Since no exact solutions of Ein- 
stein's equations are available which describe a binary in- 
spiral, we have to resort to some approximation scheme. 
The templates used by LIGO are calculated using a post 
Newtonian (pN) expansion in the orbital velocity of the 
binary However, since the pN expansion converges 

slowly U , if at all, there is a need to ascertain the quality 
of these templates. 

How does one assess the quality of these templates? 
In a first paper J?]] we compared the 2-pN templates to 
waveforms calculated from black hole perturbation the- 
ory J8|,[| . It was found that the 2-pN templates are suffi- 
ciently accurate for the detection of a binary signal, but 
not for parameter estimation. However, this first pa- 
per did not elaborate on the validity of the analysis pre- 
sented. In this publication we address this question, and 
also present results for a wider selection of detector noise 



curves. 

The paper is organized as follows. In section || we in- 
troduce the notation and briefly present the formalism 
used in the remainder of the paper. Specifically we in- 
troduce the fitting factor jnj] which is a measure of the 
quality of a template set. In section III we address some 
of the problems encountered when calculating the fitting 
factor using black hole perturbation theory. Finally, in 
section |y| we present our results, and finally in section 
^ we summarize our conclusions. 



II. ASSESSING THE QUALITY OF A TEMPLATE 

SET 

A. The fitting factor 

Let us consider the output of a gravitational wave de- 
tector of the form 



s(t) = h(t,p) + n(t). 



(1) 



Here nit) denotes the detector noise and h(t, jl) is the 
(possibly absent) signal, which depends on a set of pa- 
rameters jl. We define the signal to noise ratio (SNR) 
]Tl| associated with the template U by 



Pi = 



(s\U) 



rms(n\ti) ' 

where the inner product (-|-) is defined by 
df 



(2) 



(s\t) := 2 



S n (f) 



Hf)i*(f) + s*(f)Kf))- (3) 



As usual hats denote the Fourier transform s(f) = 
J e~ 2vi f t s(t)dt and an asterix indicates complex conjuga- 
tion. The function S n (f) is the noise's one sided spectral 
density. 

In this work we use an a naly tic fit to the noise curves 
of the various detectors |l^,[l3| : 
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for f > f r and S n (f) — oo otherwise. Table | lists the 
parameters for a few different detectors. 



Detector 


So 




Si [Hz" 1 ] 


/m[Hz] 


A [Hz] 


A [Hz] 


Caltech 40m 


0.5 x 10" 


-4-2 


10.0 


610.0 


183.0 


140.0 


initial LIGO 


0.5 x 10 


-46 


4.0 


200.0 


0.0 


40.0 


advanced LIGO 


1.5 x 10 


-49 


4.0 


70.0 


0.0 


10.0 



TABLE I. The parameters used in the analytic fit of the 
various detector noise curves. Note that the value of So is not 
important for the calculation of the fitting factor. 



For Gaussian noise 
SNR is given by 

(ft) 



LlJ the expectation value of the 



(h\U) 
y/(ti\ti)' 



(5) 



This expression is maximized for U = h, i.e. for the 
Wiener optimal filter. We can thus rewrite equation (|5|) 

as 



(Pi) 



(h\U) 



(P)r 



(G) 



The quantity 



A 



(h\U) 



V(ti\U)(h\h) 



is called the ambiguity function. It gives the loss of SNR 
due to an inaccurate template U. Its maximum FF 
max, Ai is the fitting factor JTo| . It gives the loss of SNR 
due to a non optimal set of templates. 

Obviously to calculate the fitting factor we need to 
know the signal h we are looking for. In this paper we use 
a signal calculated from black hole perturbation theory as 
a substitute for the "real" signal h || . Black hole pertur- 
bation theory is a priori only applicable in the small mass 
ratio limit. However, we believe that black hole pertur- 
bation theory can give useful information for larger mass 
ratios. We elaborate on this statement below. 

Let us now discuss the signals and templates in some 
detail. 



refer to a Cartesian system centered at the black hole 
with the z-axis orthogonal to the orbital plane. The or- 
bital velocity v is related to the angular velocity by 
the relation v 3 = Mil, where M = mi + mi is the total 
mass. 

This waveform describes a particle moving in a fixed 
circular orbit around a Schwarzschild black hole. We can 
easily generalize this expression to an adiabatic inspiral 
(see Ref. for details). To simplify expression (Q) we 
choose ip = (since the orbit is circular this represents 
no loss of generality) and a detector which singles out the 
+ polarization. Finally, we calculate the Fourier trans- 
form by means of a stationary phase approximation. We 
obtain 



l—m 



Hf) °c ^ ^ B Wl 

1=2 m=2 



(v)e im ^S lm (ti). 



(8) 



The functions Bi m (y) and ijj(v) are calculated numeri- 
cally (see Ref. @ for details) and Sj m (i?) = _ 2 ^ m (i?,0) + 
(— iy_2Yi- m ($, 0). In this expression the gravitational 
wave frequency is related to the orbital velocity by / = 
mv 3 /2irM . Black hole perturbation theory predicts that 
the innermost stable circular orbit (ISCO) is located at 
«isco = l/\/6j after which h(v) = 0. In principle, the 
sum over the index I extends to infinity. We have numer- 
ically verified that in practice we can truncate the series 
at I = 3. 



C. The templates 

As was mentioned previously, the templates to be used 
for LIGO data analysis are the restricted 2-pN waveforms 
||. (The term restricted means that we keep the phase 
accurate to second pN order but use only the leading 
order amplitude.) Using the stationary phase approxi- 
mation gives 



U :=i(/,0)oc/-*e 



«/>(/) 



(9) 



where ip(f) — 27r/t c — <f> c + </>(/) and 



B. The reference signal 

Black hole perturbation theory || gives the two inde- 
pendent components of the gravitational wave field as 

h+(t)-ih x (t)ocJ2 J2 A lm {v)e mint ^Y lm {d, V ). (7) 

1=2 rn=-l 

Here the functions Ai m (v) are the mode amplitudes 
which are calculated numerically by solving the Teukol- 
sky equation. The functions _2i^m(i?, (p) are spin- 
weighted spherical harmonics [M; the angles d and ip 



^(/) = ^gM4/r 5/3 
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/ 3058673 5429 



n 



743 
336 
617 



— 71 I x z — \Qttx 3 
4 ' 



V 1016064 1008 ' 144 



(10) 



with a; 3 = nri~ 3 / 5 Mj '. Here, 9 = {(f> c ,t c ,M,r)} are 
the template parameters, with <p c the phase at coales- 
cence (formally / = oo), t c the time at coalescence, 
M. = (mim2) ! M _ 5 the chirp mass, and rj = mim^/M 2 
the mass-ratio parameter. Note that rj is maximal for for 
mi = m2 with a value of 0.25. 
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III. VALIDITY OF THE TEST MASS 
APPROXIMATION 

Black hole perturbation theory is valid only for small 
values of -q, say rj < ?7 max = 0.1. This poses a problem, 
since most systems we are interested in have a mass ratio 
of rj > 0.2. Furthermore, for a fixed chirp mass, binaries 
with small mass ratios are more relativistic and we expect 
the pN templates to perform poorly. 

We can estimate of the minimal value that rj should 
take for the pN approximation to be valid. We expect pN 
corrections in the phase ( |l0| ) to be small if x < xq — 0.25. 
Then, for fixed chirpmass M (i.e. for a fixed number of 
wave cycles), we find 



5 



at a given frequency /. For example, for Ai = IMq, 
/ = 200 Hz (the peak sensitivity of the initial LIGO 
detector) and xq = 0.25 we get t] m i n = 0.06. Thus for 
r\ £ [0.06, 0.1] we expect out analysis to give valid results. 



0.75 1 — 

0.015 



FIG. 1. shows the fitting factor as a function of ?? ac tuai 
(i.e. the signal's rj) for constant values of the chirp mass 
M for the 40 meter noise curve. The three systems shown 
in this plot correspond to chirp masses Ma = O.87M0, 
Mb = 1.22M Q and M c = 3.92M . The primed labels cor- 
respond to signals which include all modes I — 2, 3, \m\ < I, 
whereas the unprimed labels correspond to signals only con- 
sisting of the (dominant) I = m = 2 mode. The arrows in- 
dicate where the respective curves leave the targeted mass 
range 0.8M Q < m < 3.2JW© of the detector. 



the two types of templates. This strongly suggests that 
the non-zero corrections in ( |i"o| ) play no essential role, 
and that our analysis stays valid even for rj > rj max . 

Since the binary system becomes less relativistic with 
increasing rj and constant chirp mass we expect the fit- 
ting factor to increase as well. This provides us with an 
additional consistency check. This behaviour is indeed 
observed (see Figure [l]). 



IV. RESULTS 

Calculating the fitting factor is in principle straightfor- 
ward. However evaluating the inner product (^) numer- 
ically can be tricky. Two methods were used to numer- 
ically calculate the fitting factor. In the first we calcu- 
lated the integral in (||) using a Romberg integration jl5| . 
The second method uses the fact that h(f) oc exp27rf c / 
and thus the overlap integral (J3J) can be evaluated using 
FFTs. Signal and templates in the latter case were cal- 
culated using the GRASP library Jltj} ]. For all relevant 
cases the two methods give the same answer within a few 
percent. As it turns out, the second method is a lot faster 
despite the many more "t c " values that are calculated. 

The ambiguity function depends a priory on the four 
parameters <j) c , t Cl M. and rj. Maximization over the phase 
4> c is trivial. If FFTs are used maximization over t c is also 
s trivial. This leaves us with a reduced ambiguity func- 
tion to maximize, depending on three or two variables, 
respectively. In the first case maximization is difficult, 
because A(t c , M, 77) possesses many local maxima. These 
seem to mostly disappear after maximization over t c , so 
maximizing the ambiguity function using FFTs is much 
easier. 

In a first run, reported in 0, we calculated fitting fac- 
tors for the restricted 2 pN templates assuming the ad- 
vanced LIGO noise curve. We found that the templates 
are generally adequate for detection, but not for param- 
eter estimation. We now have extended the analysis to 
also cover the Caltech 40 meter prototype and the initial 
LIGO noise curves. Our results are summarized in Ta- 
bic H for a few selected cases. For each detector there is 
at least one system with a fitting factor larger than 90% 
which is deemed acceptable by Apostolatos pel . 



However we still cannot say anything about systems 
with mass ratios larger than rj > 0.1. Furthermore, for 
some noise curves ?? m in > f?max- For example, at the peak 
sensitivity of the 40-meter noise curve (/ = 610 Hz) we 
get 7? min > 0.25. 

To estimate how important the corrections to the fit- 
ting factor due to non-zero rj are, we can set r\ = in 
the coefficients of x in the phase ([l0|). This essentially 
gives us a 2-pN approximation to the perturbation theory 
signal. Somewhat surprisingly we find, within numerical 
accuracy, no difference in the fitting factors calculated for 





Advanced LIGO 


initial LIGO 


40-meter 


Templates 


Newtonian 


2 pN 


2 pN 


2 pN 


1.4 - 1.4M 


97.2% 


93.0% 


95.8% 


94.0% 


0.5 - 5.OM 


51.6% 


95.2% 


89.7% 


89.0% 


1.4 - 10.0Af Q 


55.8% 


91.4% 


88.4% 


< 85.0% 


10.0 - 10.0M Q 


70.1% 


91.7% 


86.3% 




4.0 - 3O.OM 


61.3% 


86.6% 


67.9% 





TABLE II. The fitting factors FF for selected binary sys- 
tems and for various noise curves. 
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A more complete picture is given by Figure || which 
shows the contours of FF = 90% in the mi — m 2 plane 
for the 40 meter, initial and advanced LIGO detectors. 
Since the frequency band for the advanced LIGO starts 
at 10Hz we have to follow a large number of wave cycles, 
especially for low mass systems. This makes the evalu- 
ation of the fitting factor difficult in terms of available 
computer power. 
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FIG. 3. shows plots of the fitting factor FF as a func- 
tion of r] with constant chirp mass M for the same three 
systems. The curves shown were calculated using test mass 
signals consisting of the I = m = 2 mode only. The dashed 
curves correspond to templates with r\ = 0. The solid lines 
represent proper 2 pN templates. There is practically no dif- 
ference between the two cases, indicating that perturbation 
theory gives valid results even for relatively large mass ratios. 



FIG. 2. shows the contours of FF = 90% in the mi — mi 
plane for the Caltech 40 meter prototype, the initial and ad- 
vanced LIGO detectors. Units are in solar masses. The dotted 
curve to the right denotes the line where the frequency at the 
inner most stable orbit (ISCO), as predicted by perturbation 
theory, becomes less than 100Hz. While there are still tem- 
plates with FF > 90% in this region, the total SNR for these 
systems is very low since they spend very little time in the 
advanced LIGO frequency band. 

In a next step we calculated the fitting factor as func- 
tion of rj for fixed chirp mass M.. (These are the mass 
parameters for the signal and they should not be confused 
with the mass parameters in the ambiguity function A.) 
As expected, the fitting factor drops considerably as ij 
decreases. Figure [l] shows the curves FF(rj) for three 
different values of the chirp mass for the 40 meter noise 
curve. Note that the end points of these curves lie well 
outside the mass regime considered of interest for this 
detector. Similar results were obtained for the advanced 
and initial LIGO detectors. 

Finally, we repeated this calculation with rj = in 
the coefficients of the phase (10). The values for the 



fitting factor were found to agree with the previous values 
within numerical accuracy. Figure [| compares the two 
cases for the 40m noise curve for the I = m = 2 signal. 
We take this as an indication that perturbation theory 
waveforms are adequate to test the performance of a set 
of templates, even when rj is not small. 



V. CONCLUSIONS 

The results of this analysis suggest that the 2 pN tem- 
plates are adequate for the detection of gravitational ra- 



diation emitted from inspiraling compact binary systems, 
when the masses lie in the region shown in Figure ||. 
We emphasize that we do not claim that the 2 pN tem- 
plates are all we need to know for a LIGO-type search 
for compact binaries. Clearly, more accurate templates 
are needed to detect the larger mass systems as well as 
for parameter estimation JtJ. But we do believe that the 
2-pN templates give us a good chance to detect gravita- 
tional wave from, at least, some systems. 

While it is not a priory clear that black hole pertur- 
bation theory makes a good tool to ascertain the quality 
of 2-pN templates for large mass ratio systems, we have 
given arguments to support our results. At the very least, 
we feel our analysis gives a qualitatively correct picture. 



ACKNOWLEDGMENTS 

The authour would like to thank Eric Poisson for his 
stimulating comments. This work was supported by the 
Natural Sciences and Engineering Research Council of 
Canada. 



[1] K. S. Thorne, Gravitational radiation, in Three hundred 
years of gravitation, edited by S. W. Hawking and W. Is- 
rael, chap. 9, pp. 330-458, Cambridge University Press, 
Cambridge, 1987. 

[2] L. A. Wainstein and V. D. Zubakov, (Prentice-Hall, Lon- 
don, 1962). 



4 



L. Blanchet, B. R. Iyer, C. M. Will, and A. G. Wise man, 
Class. Quant. Grav. 13, 575 (1996), |gr-qc/9602024 . 
C. M. Will and A. G. Wiseman, Phys. Rev. D54, 4813 



(1996), kr-qc/9608012 . 

L. Blanchet, T. D amour, and B. R. Iyer, Phys. Rev. D51 
5360 (1995), |gr-qc/9501029. 



E. Poisson, Phys. Rev. D52, 5719 (1995), |gr-qc/950503C . 
S. Droz and E . Poisson, Phys. Rev. D56, 4449 (1997), 
gr-qc/9705034 . 

E. Poisson, Phys. Rev. D47, 1497 (1993). 

S. A. Teukolsky, APJ 185, 635 (1973). 

T. A. Apostolatos, Phys. Rev. D54, 2421 (1996). 

C. Cutl er and E. E. Flanagan, Phys. Rev. D49, 2658 



(1994), kr-qc/9402014 



A. Abramovici et al, Science 256, 325 (1992). 

E. E. F lanagan and S. A. Hughes, Phys. Rev. D57, 4535 
(1998), kr-qc/970103~ . 

J. Goldberg, A. MacFarlane, E. Newman, F. Rohrlich, 
and E. Sudarshan, J. Math. Phys. 8, 2155 (1967). 
W. A. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. 
Flannery, Numerical Recipes: The Art of Scientific Com- 
puting (Cambridge University Press, Cambridge, 1995). 

B. Allen, GRASP: a data analysis package for grav- 
itational wave detection (LIGO Collaboration, 1997), 



http:/ /www. ligo.caltech.edu/LIGO_web/Collaboration/manual.pdf 



5 



